Constraints on neutron star radii based on chiral effective field theory interactions 
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We show that microscopic calculations based on chiral effective field theory interactions constrain 
the properties of neutron-rich matter below nuclear densities to a much higher degree than is reflected 
in commonly used equations of state. Combined with observed neutron star masses, our results lead 
to a radius R = 9.7 — 13.9 km for a 1.4 Mg star, where the theoretical range is due, in about equal 
amounts, to uncertainties in many-body forces and to the extrapolation to high densities. 
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Introduction. - With the advances in observational ca- 
pabilities, it is becoming possible to obtain direct evi- 
dence for the size of neutron stars [l| . Sources of infor- 
mation include measurements of optical radiation from 
nearby isolated neutron stars whose distances are known 
from parallax determinations, and observations of ther- 
monuclear flares on the surfaces of neutron stars [2, [3j ■ 
Neutron star seismology [4|, pulse profiles in X-ray pul- 
sars [5| and moment of inertia measurements [6| are ad- 
ditional sources. In the near future, one further ex- 
pects gravitational wave signals from collapsing stars 
and merging binary neutron stars to provide information 
about the equation of state of dense matter [7j . 

In nuclear physics, recent developments of effective 
field theory (EFT) and the renormalization group (RG) 
for nuclear forces enable controlled calculations of matter 
at nuclear densities [8j4l0|. In this framework, it is pos- 
sible to estimate the theoretical uncertainties due to ne- 
glected many-body forces and from an incomplete many- 
body calculation. In this Letter, we show that such cal- 
culations of the equation of state (EOS) below nuclear 
densities constrain the properties of dense matter, and 
the radii of typical neutron stars, to a much higher de- 
gree than is reflected in current neutron star modeling. 

Neutron matter below nuclear densities.- Our studies 
are based on microscopic calculations of neutron matter 
using chiral low-momentum two-nucleon (NN) and three- 
nucleon (3N) interactions [ll| . The largest uncertainty 
for the neutron matter energy arises from the couplings 
C3 and (to a smaller extent) c\ that determine the lead- 
ing two-pion-exchange three-body forces among neutrons 
in chiral EFT. We improve the range of C3 compared to 
Ref. [ill by taking C3 values from extractions based on 
the same couplings in the subleading two-pion-exchange 



rj-UJ. In ad- 



NN interaction, c 3 = -(3.2 - 4.8) GeV 
dition, we include a shift 5c^ = 1.0 GeV -1 to take into 
account contributions at the next order for 3N forces [8j . 
In the following we therefore use C3 = —(2.2 — 4.8) GeV - 



and C\ = —(0.7—1.4) GeV with a smooth 
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ulator. The neutron matter calculations details are as in 
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Ref. 1 1] . which suggested that the energy is perturbative 
at nuclear densities. Using only NN interactions, we ob- 
tain a neutron matter energy per particle E™ (po)/N — 
10.4 MeV at saturation density p — 2.7 x 10 14 gcm~ 3 . 
(We define the density p as the nucleon mass times the 
baryon density.) The inclusion of 3N forces leads to 
E n (p )/N = 16.3 ± 2.2 MeV, dominated by the repul- 
sive C3 contribution and the associated uncertainty. The 
3N contribution of « 6 MeV is to be compared to the 
NN potential energy ra —26 MeV (the kinetic energy is 
3/5£f(/9o) ~ 36 MeV), and also the 3N uncertainty of 
« 2 MeV is consistent with the contributions of higher- 
order 3N forces (given an approximate expansion param- 
eter in chiral EFT of ~ 1/3 at these momenta). Other 
microscopic calculations lie within our theoretical uncer- 
tainties [llj, and at low densities p ~ po/10, the results 
are also consistent with calculations for resonant Fermi 
gases including effective range contributions [15j . 

Neutron star matter. - We extend our results to matter 
in beta equilibrium using the parametrization: 
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(1) 

where E n (p)/N is given by our neutron matter results, 
to is the nucleon mass, and x the proton fraction. The 
last terms in Eq. (JTJ) incorporate the contributions from 
protons through the symmetry energy S 2 (p) and from 
electrons [161 ] . The proton fraction in beta equilibrium is 
given approximately by x(p) « [4S , 2(p)/(fi.c)] 3 /(37r 2 p/77i), 
and for S 2 (po) ~ 30 MeV, x(p ) ~ 5%. The energy dif- 
ference between neutron matter and matter in beta equi- 
librium is —1.1 MeV, and this will have only a minor 
impact on our final results. We extract S 2 {p) for nuclear 
densities using empirical saturation properties, 
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FIG. 1: Pressure of neutron star matter based on chiral low-momentum interactions for densities p < pi (corresponding 
to a neutron density pi, n = l.lpo). The band estimates the theoretical uncertainties from many-body forces and from an 
incomplete many-body calculation. At low densities, the results are compared to a standard crust EOS [12J], where the right 
panel demonstrates the importance of 3N forces. The extension to higher densities using piecewise polytropes (as explained in 
the text) is illustrated schematically in the left panel. 



with binding energy ay = 16MeV and incompressibility 
K = 230 MeV (which are within theoretical uncertainties 
of the nuclear matter calculations of Ref. [10(). To in- 
clude the symmetry energy in Eq. ([TJ , we use the Ansatz 
^(p) = <S l 2(p/po) 7 and fit S2 = ^(po) an d 7 to our neu- 
tron matter results. The fit has a relative uncertainty of 
< 5% for densities p /8 < p < pi — 3.0 x 10 14 gcm~ 3 
(pi corresponds to a neutron density p\^ n = l.lpo). We 
obtain the following symmetry energy parameters and 
proton fractions: 
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The resulting pressure of neutron star matter is shown 
in Fig. [TJfor densities p < pi, where the band is domi- 
nated by the uncertainty in C3. The comparison of these 
parameter-free calculations to a standard crust EOS 17] 
shows good agreement to low densities p > po/10 within 
the theoretical uncertainties. In addition, the right panel 
of Fig. [Tj demonstrates the importance of 3N forces. The 
pressure obtained from low-momentum NN interactions 
only, based on the RG-evolved N 3 LO potentials of Entem 
and Machleidt (EM) [13] or Epelbaum et al. (EGM) [lj|, 
differ significantly from the crust EOS at pa/2. 

Neutron stars. - The structure of non-rotating neutron 
stars without magnetic fields is determined by solving 
the Tolman-Oppenheimer-Volkov (TOV) equations. Be- 
cause the central densities reach values higher than p\, 
we need to extend the uncertainty band for the pressure 
of neutron star matter beyond p\ . To this end, we intro- 
duce a transition density p i2 that separates two higher- 
density regions, and describe the pressure by piecewise 
polytropes, P(p) — n\p Tl for p\ < p < p±2, and 



P(p) = K2p T ' 2 for p > pi2, where Ki,2 are determined by 
continuity of the pressure. Ref. [18] has shown that such 
an EOS with 1.5 < Ti.2 < 4.0 and transition densities 
P12 ~ (2.0 — 3.5)po can mimic a large set of neutron star 
matter EOS. We therefore extend the pressure of neu- 
tron star matter based on chiral EFT in this way, with 
1.5 < I\2 < 4.5 and 1.5 < pu/po < 4.5, as illustrated 
in Fig. [TJ The possibility of a phase transition at higher 
densities is implicitly taken into account if one regards 
the 1^2 values as averages over some density range. 

We solve the TOV equations for the limits of the pres- 
sure band below nuclear densities continued to higher 
densities by the piecewise polytropes. The range of r 12 
and P12 can be constrained further, first, by causal- 
ity, which limits the sound speed to the speed of light 
and, second, by the requirement that the EOS support 
a 1.65M0 star [19J] . The resulting allowed range of poly- 
tropes is shown by the light blue band at higher density 
in Fig. [5] [27J. The comparison with a representative set 
of EOS used in the literature [16j demonstrates that the 
pressure based on chiral EFT interactions (the darker 
blue band) sets the scale for the allowed higher-density 
extensions and is therefore extremely important. It also 
significantly reduces the spread of the pressure at nu- 
clear densities from a factor 6 at p\ in current neutron 
star modeling to a factor 1.5. 

Results.- In Fig. [3] we show the neutron star M-R 
curves obtained from the allowed EOS range. The blue 
region corresponds to the blue band for the pressure in 
Figs. Q] and [2] At the limits of this region, the pressure 
of neutron star matter is continued as piecewise poly- 
tropes, and all curves end when causality is violated. 
Should causality be violated before the maximum mass 
(at dM/dR = 0) is reached, one could continue the M- 
R curves by enforcing causality. This would lead to a 
somewhat larger maximum mass, but would not affect 
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FIG. 2: Comparison of the EOS based on Fig.Q]to a represen- 
tative set of EOS used in the literature [16j]. The blue band 
corresponds to the band in Fig.[T]and the lighter region covers 
the range of polytropes allowed (see text for discussion). 



the masses and radii of neutron stars with lower central 
densities. We observe from the transition density points 
P12 in Fig. [3] that the range of Ti dominates the uncer- 
tainty of the general extension to high densities. Smaller 
values of Ti are excluded because the associated EOS is 
not able to support a 1.65Af Q star. The larger allowed 
values of the polytropic indices lead to very low central 
densities p ~ (2.0 — 2.5)po- 

We find that the pressure at nuclear densities and be- 
low sets the scale for the M-R results. The blue region 
in Fig. [3] ends almost at the central value of the radius 
results. For a 1.4M Q star, the radius is constrained to 
R = 9.3 — 13.5 km, as indicated by the vertical band. 
While going from neutron matter to beta equilibrium 
can reduce individual results for an IAMq star by up to 
0.4 km, the overall result is very similar for pure neutron 
matter with R = 9.3 — 13.3 km. Furthermore, if a 2.OM 
star were to be observed, this would reduce the allowed 
range to R = 10.5 - 13.3 km. As for the EOS in Fig. H 
the presented radius constraint significantly reduces the 
spread of viable neutron star models, e.g., it is difficult 
to see how one can obtain R s=s 15 km as is the case for 
the Shen EOS [20(. Finally, our results are more rigorous 
than an estimate based on the empirical PR~^ correla- 
tion [16j |. which for the values of the pressure we find, 
P{po) = 1.4-2.1MeVfm~ 3 , implies R = 9.4- 11.9km. 

When chiral 3N forces are neglected, the neutron star 
radius is significantly smaller, with i? NN = 8.8 — 11.0 km 
as shown in Fig. |4] based on low- momentum NN interac- 
tions only. This demonstrates that the theoretical error 
for the radius of a 1.4Mq star is due, in about equal 
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FIG. 3: Neutron star M-R results for the EOS based on 
Fig. [l] The thick (thin) lines, corresponding to the left (right) 
branch, start from the low pressure limit ci = —0.7 GeV - , 
c 3 = -2.2 GeV -1 (high pressure limit ci = -1.4 GeV -1 , 
d — —4.8 GeV -1 ). The blue region corresponds to the band 
below nuclear densities in Figs. [T] and [21 The different piece- 
wise polytropes can be identified from the colors/lines indi- 
cating Ti/r2 and from the points denoting pi2. The vertical 
band gives the radius constraint for a 1.4Mg star. 



amounts, to the uncertainties in 3N forces and to the 
extension to higher densities dominated by T\. 

Effect of the crust.- In our calculations, the difference 
between the neutron and proton masses was neglected 
and the phases were assumed to be spatially uniform. 
In this approximation, matter at low density consists 
only of neutrons. The impact of using a more realis- 
tic EOS at low densities can be investigated by observ- 
ing that the surface gravity of the star is approximately 
constant in the outer layers. By integrating the equa- 
tion of hydrostatic equilibrium from the surface of the 
star up to a crust density p c , one finds that the mass 
between the density p c and the surface is proportional 
to the pressure at p c 



21] , Thus the stellar mass is to a 



good approximation unaffected by changes in the EOS 
at p < p c . To determine how changes in the low-density 
EOS affect the radius, we note that the thickness of the 
crust {p < p c ) is AR — [p(p c ) ~ Hs\/[mg(l + z)], where 
g = GM(1 + z)/R 2 is the surface gravity, with surface 
redshift 1 + 2 = [l-2GM/{Rc 2 )}- 1 / 2 p\. Here p s is the 
(neutron) chemical potential at the surface of the star, 
where the pressure is zero. For the calculations in this 
paper, p s — mc 2 , while for realistic EOS of cold catal- 
ysed matter it includes the binding energy per particle 
of solid iron, « 8MeV. Thus use of a more realistic 
EOS at low densities will increase the radius of the star 
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FIG. 4: Same as Fig. [3] but for NN-only interactions (thick 
lines are based on EM [12| and thin lines on EGM U3 



z)]. This increases the radius for a 
— 0.5 km, leading to our final result 



by 8MeV/[mg(l + 
1AM & star by 0.2 
# = 9.7- 13.9 km. 

Other implications.- The relatively weak density de- 
pendence of the nuclear symmetry energy also makes 
predictions for the neutron skin of 208 Pb. The symme- 
try energy of a nucleus in the liquid droplet model con- 
sists largely of bulk and surface contributions, the lat- 
ter being determined by an integration of 52 (p) through 
the nucleus [23j. Assuming a quadratic density de- 
pendence for the energy of symmetric nuclear matter 
and our results for 52 = S2(p/po) J , one finds the ra- 
tio of the surface and bulk symmetry parameters to be 
S s /S2 ~ 1-85 ± 0.25. This leads to a neutron skin thick- 
ness 5R = 0.16 — 0.2 fm for 208 Pb, while the correlation 
with the slope of the neutron matter energy 24| gives 
0.14- 0.19 fm. Therefore we predict SR = 0.14 - 0.2 fm, 
which can be tested in the parity- violating electron scat- 
tering experiment (PREX) [25|. Finally, we note that in 
a complementary approach [26|, the EOS of high-density 
matter is constrained from perturbative QCD calcula- 
tions, and that our results are very consistent with the 
astrophysical estimates of Ref . [3J . 

In this Letter, we have demonstrated that microscopic 
calculations based on chiral EFT and many-body theory 
constrain the pressure of matter at nuclear densities to 
within ±25%. This should be taken into account in mod- 
eling stellar collapse, black hole formation, and neutron 
stars. Even allowing for uncertainties in the low-energy 
theory and the extrapolation to higher densities, we find 
that the radius of a neutron star depends only weakly its 
mass and for a 1.4M Q star is rather well constrained. 
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